Searching for Dark Matter in the CMB: 
A Compact Parameterization of Energy Injection from New Physics 



Douglas P. Finkbeiner, 1 ' 2 ^ Silvia Galli, 3 ' 4 <0 Tongyan Lin, 2 '! and Tracy R. Slatyer 5 '! 

Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 
''Physics Department, Harvard University, Cambridge, MA 02138, USA 
Laboratoire Astroparticule et Cosmologie (APC), 
Universite Paris Diderot, 75205 Paris cedex 13, France 
Physics Department and INFN, Universita' di Roma "La Sapienza", Pie Aldo Moro 2, 00185, Rome, Italy 
J School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA 

High-precision measurements of the temperature and polarization anisotropies of the cosmic mi- 
crowave background radiation have been employed to set robust constraints on dark matter annihi- 
lation during recombination. In this work we improve and generalize these constraints to apply to 
energy deposition during the recombination era with arbitrary redshift dependence. Our approach 
also provides more rigorous and model-independent bounds on dark matter annihilation and decay 
scenarios. We employ principal component analysis to identify a basis of weighting functions for the 
energy deposition. The coefficients of these weighting functions parameterize any energy deposition 
model and can be constrained directly by experiment. For generic energy deposition histories that 
are currently allowed by WMAP 7 data, up to 3 principal component coefficients are measurable 
by Planck and up to 5 coefficients are measurable by an ideal cosmic variance limited experiment. 
For WIMP dark matter, our analysis demonstrates that the effect on the CMB is described well by 
a single (normalization) parameter and a "universal" redshift dependence for the energy deposition 
history. We give WMAP 7 constraints on both generic energy deposition histories and the universal 
WIMP case. 

PACS numbers: 95.35.+d,98.80.Es 



I. INTRODUCTION 

Measurements of the cosmic microwave background 
(CMB) in the past decade by experiments including 
WMAP, ACBAR and BOOMERANG [H-Q, and more 
recently SPT, QUaD and ACT have provided an 

unprecedented window onto the universe around redshift 
1000. With the advent of the Planck Surveyor 0, the 
successor experiment to WMAP, percent-level modifica- 
tions to recombination will be observable. Planck has 
already completed three sky surveys and begun a fourth, 
and cosmological data are expected to be released pub- 
licly in 2012-13. 

Accurate measurements of the CMB have the poten- 
tial to probe the physics of dark matter (DM) beyond its 
gravitational interactions. In the large class of models 
where the DM is a thermal relic, its cosmological abun- 
dance is determined by its annihilation rate in the early 
universe: the correct relic density (~ 22% of the energy 
density of the universe) is obtained for an s-wave anni- 
hilation cross section of (crv) f=s 3 x 10~ 26 cm 3 /s during 
freezeout. 

DM annihilation at this rate modifies the ionization 
history of the universe and has a potentially measurable 
effect on the CMB. During the epoch of recombination, 
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DM annihilation produces high-energy photons and elec- 
trons, which heat and ionize the hydrogen and helium 
gas as they cool. The result is an increased residual ion- 
ization fraction after recombination, giving rise to a low- 
redshift tail in the last scattering surface. The broader 
last scattering surface damps correlations between tem- 
perature fluctuations, while enhancing low-£ correlations 
between polarization fluctuations. 

The resulting constraints on the dark matter annihi- 
lation rate have been studied by several authors [cj-[T3|. 
These bounds have a notable advantage over other indi- 
rect constraints on dark matter annihilation, in that they 
are independent of the DM distribution in the present 
day, and do not suffer from uncertainties associated with 
Galactic astrophysics. They depend only on the cosmo- 
logical DM density, which is well measured; the DM mass; 
the annihilation rates to the final states; and the stan- 
dard physics of recombination. Recombination modeling, 
while not simple, involves only well-understood conven- 
tional physics, and the latest models are thought to be 
accurate at the sub-percent level required for Planck data 

Current limits from WMAP already significantly con- 
strain models of light dark matter with masses of around 
a few GeV and below, if the annihilation rate at recom- 
bination is the thermal relic cross section. Heavier DM 
is also constrained if the annihilation rate is enhanced 
at low velocities or for other reasons is much larger than 
the thermal relic cross section at recombination. Mod- 
els lying in these general categories are also of significant 
interest for their possible connection with experimental 
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anomalies. 

The DAM A/LIBRA [17] and CoGeNT QJ, HU di- 
rect detection experiments have reported excess events 
and annual modulation that may have a consistent ex- 
planation as originating from scattering of 0(5 — 10) 
GeV WIMPs (e.g. [20]). Results from the XENON10, 
XENON100 and CDMS experiments are in tension with 
this interpretation [2ll [22| , but there is ongoing debate 
on the sensitivity of these experiments to the very low- 
energy nuclear recoils in question (see e.g. [H, EH)- 

The PAMELA [H, Fermi Sl> PPB-BETS [13, 
ATIC [28] and H.E.S.S have measured electron and 
positron cosmic rays in the neighborhood of the Earth, 
and found results consistent with a new primary source 
of e ± in the 10 — 1000 GeV energy range. If the sig- 
nal is attributed to dark matter annihilation then the 
annihilation rate in the Galactic halo today must be 1- 
3 orders of magnitude above the canonical thermal relic 
value [13, HJ . This has motivated models of dark mat- 
ter with enhanced annihilation at low velocities [HI, [33[ . 
While this enhancement would not be significant during 
freezeout, it would be effective at recombination when 
the typical velocity of dark matter is v ~ 10 _8 c Q. 

With the release of data from Planck expected in the 
next two years, models falling into these categories should 
either be robustly ruled out, or give rise to a measurable 
signal d, [HJ EH • If no signal is observed, the sensitivity 
of Planck will allow us to probe regions of parameter 
space relevant for supersymmetric models, where the DM 
is a thermal relic with mass of several tens of GeV. It is 
timely to explore improvements to these constraints. 

The approach of previous studies has been to specify 
the energy deposition history (redshift dependence) and 
then calculate the effect on the ionization history and 
anisotropy spectrum using public codes such as RECFAST 
and CAMB. A single parameter describing the normaliza- 
tion of the signal is then added to the standard likelihood 
analysis using CosmoMC, and bounded by WMAP obser- 
vations. The redshift dependence has been studied in 
two cases: in the "on-the-spot" case, assuming that the 
amount of energy deposited to the gas precisely tracks 
the rate of dark matter annihilation (e.g. @, [9lll3|). or 
employing detailed energy deposition histories for spe- 
cific models (e.g. [13, 13113) • 

In the first case, model- 
independent constraints are obtained, but without a pre- 
cise way to connect the bounds to any particular model. 
The second case only precisely constrains specific models. 

While these analyses have been adequate for simple 
estimates of whether a model is strongly ruled out, eas- 
ily allowed, or on the borderline, the upcoming high- 
precision data from Planck demand a more careful model- 
independent analysis. Such an analysis can also be ap- 
plied to more general classes of energy deposition his- 
tories during and after recombination: for example, the 
energy deposited by a late-decaying particle species, de- 
cay from an excited state of the dark matter, or dark 
matter annihilation in models where the redshift depen- 
dence of the annihilation rate has an unusual form (as in 



some models of asymmetric dark matter). 

In this work we exploit the fact that the effects of 
energy deposition at different redshifts are not uncor- 
rclatcd. Any arbitrary energy deposition history can be 
decomposed into a linear combination of orthogonal ba- 
sis vectors, with orthogonal effects on the observed CMB 
power spectra (CVs). For a broad range of smooth energy 
deposition histories, the vast majority of the effect on the 
Cgs can be described by a small number of independent 
parameters, corresponding to the coefficients of the first 
few vectors in a well-chosen basis. These parameters in 
turn can be expressed as (orthogonal) weighted averages 
of the energy deposition history over redshift. 

We employ principal component analysis (PCA) to 
make this statement quantitative and derive the relevant 
weight functions, and the corresponding perturbations 
to the Ce spectra. Our approach in principle generalizes 
to all possible energy deposition histories. To investigate 
the number of observable parameters, we consider generic 
perturbations about two physically interesting fiducial 
cases. We focus primarily on the example of dark matter 
annihilation, or any other scenario where the power de- 
posited per volume scales approximately as (1 + z) e (i.e. 
as density squared) , as an energy deposition mechanism, 
but also show results for the case of dark matter decay, 
or similar scenarios where the power deposited scales as 
(l + zf. 

Our computation of the effects of energy deposition on 
the CMB anisotropies, and the approximations we use 
for estimating the significance of these effects in exper- 
imental data, arc described in EjTTJ In £11111 we present 
our principal component analysis for both "annihilation- 
like" and "decay-like" general energy deposition histo- 
ries 1 . There are significant degeneracies between energy 
deposition and perturbations to the cosmological param- 
eters, and so we marginalize over the standard cosmologi- 
cal parameters when deriving the principal components 2 . 

We then address the constraints on and detectability 
of the principal components in current and future exper- 
iments. Given a Ce spectrum observed by an experiment 
(e.g. Planck), we can measure the residual with respect 
to the best-fit standard ACDM model, and then project 
this residual onto the CV-space directions corresponding 
to the principal components. Given any model for the 
energy deposition history, we can then ask if the recon- 
structed coefficients for the various principal components 
are consistent with the model. Of course, for the later 
principal components the effect on Cg's is so small that 



1 Files containing the results of these analyses are available online 
at http://nebel.rc.fas.harvard.edu/epsilon/; see also Ap- 
pendix \C\ 

2 We test the effect of including additional cosmological parame- 
ters (running of the scalar spectral index, the number of mass- 
less neutrino species, and the primordial He fraction) and find 
no large degeneracy with energy injection, justifying our neglect 
of these additional parameters in our main analysis. 
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very little information on their coefficients can be recov- 
ered. In QYV[ we make this statement quantitative, and 
estimate the number of principal components whose co- 
efficients could be detectable in Planck and an ideal cos- 
mic variance limited (CVL) experiment, subject to con- 
straints from WMAP 7. The CVL case presents a hard 
upper limit on the number of independent parameters de- 
scribing the energy deposition history that can profitably 
be retained in the analysis. We also discuss the bias to 
the standard cosmological parameters, in the case where 
there is non-zero energy deposition that is neglected in 
the analysis; in our framework it is straightforward to 
characterize the biases to the cosmological parameters 
for an arbitrary energy deposition history. 

In f|V]we present a separate principal component anal- 
ysis for the more limited case of conventional GeV-TeV 
WIMPs annihilating to Standard Model final states. We 
demonstrate that in this case, all the effect of dark mat- 
ter annihilation can be captured by one parameter only, 
i.e. the amplitude of the first principal component. 

Finally, in j jVT) we estimate the constraints on the prin- 
cipal components obtainable with current (WMAP 7) 
and future (Planck, CVL) experiments with a full likeli- 
hood analysis using the CosmoMC code. We employ here 
the principal components obtained with the Fisher ma- 
trix analysis - which assumes that the effect on the CMB 
scales linearly with the energy deposition. We illustrate 
the range of validity of this assumption for the different 
experiments considered. We check that the constraints 
previously obtained with our Fisher matrix analysis - 
which assumes Gaussian likelihood functions - are com- 
patible with the ones obtained with the CosmoMC analysis. 
We check that the constraints on a given energy depo- 
sition history can be reconstructed from the constraints 
on the principal components. 

Appendix [A] considers the effects on the analysis of 
changing various assumptions and conventions, includ- 
ing the effect of additional cosmological parameters and 
using different codes to calculate the ionization histories. 
We find that the only such choice that non-negligibly 
modifies the early (detectable) principal components is 
the treatment of Lyman-a photons, although the inclu- 
sion of additional cosmological parameters can change 
the constraints at the ~ 10% level. Appendix |B] dis- 
cusses marginalization over the cosmological parameters. 
Appendix [C] describes the results from this analysis that 
we have made available online. 



II. THE EFFECT OF ENERGY INJECTION 

We begin by considering DM annihilation-like or 
decay-like energy deposition histories. The energy in- 
jection from these sources scales respectively as density 
squared and density, so these cases cover the generic sce- 
narios where energy is injected by two-body or one-body 
processes. It is convenient to express the energy injection 
as a slowly varying function of z that depends on the 
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FIG. 1: Rate of Hydrogen ionization from energy deposition, 
relative to the number density of ionized Hydrogen (n° on ) 
when there is no energy deposition. The lines shown are the 
cases of constant p an n and pdec, corresponding to on-the-spot 
energy deposition from dark matter annihilation and dark 
matter decay, respectively. 



source of the energy injection (e.g. the WIMP model) 
and a term containing cosmological parameters. We pa- 
rameterize the energy deposition histories, respectively, 
as, 



dE 
dt dV 

dE 
dt dV 



= p am (z)c 2 £ll M p 2 c (l + zf , 

= Pdcc(z)c 2 D,BMPc(l + z) 3 , (1) 



dec 



where p ann (z) (or pd oc (z)) contains all of the informa- 
tion about the source of energy injection and the effi- 
ciency with which that energy ionizes the gas. We gener- 
ically refer to p ann and pdec as the "energy deposition 
yield." For consistency with [l4[, we express p ann (z) in 
units of cm 3 /s/GeV, while the units of Pdec(^) are s _1 . 
If the energy injection is due to DM annihilation, p an n = 
f(z)((Tv)/mr>M 1141 • where f(z) is an 0(1) dimensionless 
efficiency factor [lfj ; if the energy injection is due to DM 
decay, Pdcc(z) = /(z)/r, where r is the decay lifetime. 
Other authors have written p ann in units of m 3 /s/kg 
or parameterized the energy deposition in eV/s/baryon 
[1, 03 EH • For calibration, the energy deposition from a 
100 GeV thermal relic WIMP with f(z) = 1 corresponds 
to p ann « 3 x 10- 28 cm 3 /s/GeV « 1.7 x 10- 7 m 3 /s/kg, 
or an energy deposition of 2.1 x 10~ 24 eV/s/H, assum- 
ing the WMAP 7 best-fit cosmology. Throughout this 
work, we employ the cosmological parameters from [34j ] 



as a baseline: explicitly, ojb = 2.258 x 10 



0.1109, 



A s (fc=0.002 Mpc- 1 ) = 2.43 x 10~ 9 , n s = 0.963, r = 
0.088, H = 71.0 km/s/Mpc. 

Energy deposition during recombination primarily af- 
fects the CMB through additional ionizations; the mod- 
ified ionization history then leads to an increased width 
for the surface of last scattering. An alternate approach 
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to studying generic energy deposition histories might be 
to study generic ionization histories [35j j , since the former 
can be directly mapped to the latter. We frame the prob- 
lem in terms of energy deposition histories because this 
provides direct constraints on physical energy injection 
models. 

Suppose we are interested primarily in a class of en- 
ergy deposition histories for which the energy deposition 
yield p(z) (that is, p snn (z) or p& ec (z), as appropriate) is 
not very rapidly varying. Then we can discretize p(z) as 
a sum over a basis of N (5-function-likc energy deposi- 
tion histories, p(z) — Yli=i a i@i( z )- The basis functions 
Gi(z) are (by default) Gaussians with a = Az/4, cen- 
tered on Zi (i = 1.--/V), where Az is the spacing between 
the Zi. They are normalized such that J dzGi(z) = Az. 
For example, in the limit of large N an energy deposition 
history with constant p(z) — po corresponds to aj = po 
for all i. 

If the energy deposition is small enough, the ef- 
fect on the CMB anisotropy power spectrum is lin- 
ear in the energy depositions at different redshifts 
[6Q(p(z)) = SQiEti^Giiz)) = Y£ =1 8C l {a i G i (z))], 
and in the amount of energy deposition at any redshift 
[<5CHEL onGi(z)) = ai5Ci{Gi(z))]. Then the ef- 
feet of an arbitrary energy deposition history can be de- 
termined simply from studying the basis functions Gi{z). 
We will assume linearity throughout this work, and jus- 
tify that assumption in £ ]II Dl 

Of course, given any annihilation-like energy deposi- 
tion history, it can be rewritten in decay-like form with 
a strongly redshift-dependent p{z), and vice versa. The 
basis of Gi (z) functions can describe any energy deposi- 
tion history, at least in the large- N limit. However, the 
very different "underlying" redshift dependence in the 
two cases, and the uncertainties associated with the an- 
nihilation rate at low redshift (due to the onset of struc- 
ture formation), motivate us to study different redshift 
ranges in the two cases. 

For each Gi, we can compute the effect on the 
ionization history and the anisotropy spectrum in 
the limit of small energy deposition. We determine 
dCf T /da l: dC EE /da l: dCj^/da, Vi,L In our default 
analysis we employ the CosmoRec and CAMB codes, with 
the prescription for including the extra energy deposition 
laid out in [j| [36[ . If there are N basis functions and we 
take rig spherical harmonics into account, this yields an 
ri£ x N transfer matrix T whose (£, i)th element is, 



dCj 
da,. 



dCV dCf E dCj E 



da; 



dai 



da,. 



(2) 



In this work we focus primarily on annihilation-like en- 
ergy deposition histories, for which we restrict ourselves 
to the 80 < z < 1300 range; as a default, we will take 
50 redshift bins covering this range. At higher redshifts 
the universe is ionized and so the effect of energy deposi- 
tion on the ionization history is negligible, while at lower 
redshifts the DM number density becomes so small that 



the energy injected from annihilation is insignificant, as 
shown in Figure [T] This in turn justifies our neglect of 
DM structure formation: while for z < 100, DM clumps 
start to form and the annihilation rate no longer tracks 
the square of the average relic density, the energy in- 
jection is already sufficiently suppressed that the signal 
remains negligible. 

For DM decay, the signal is not nearly so suppressed at 
low redshifts, and so we consider the redshift range 10 < 
z < 1300. With this expanded redshift range, we switch 
from linear to log binning, with 90 bins covering this 
redshift range 3 ; we take the basis functions Gi(ln(l + z)) 
to be Gaussians in ln(l + z), normalized so that their 
integral with respect to dln(l + z) is given by the spacing 
between the log bins Aln(l + z). With these choices, 
again an energy deposition history with constant p(z) = 
Po corresponds (in the large- TV limit) to po J2i=i G%- 

We again ignore structure formation in the decay-like 
case, where the total power injected depends only on the 
average density. The universe is rather transparent to the 
products of DM decay and annihilation at these redshifts, 
so even a very spatially non-uniform distribution of en- 
ergy injection would not be expected to cause ionization 
or temperature hot-spots (at least for particles injected 
at weak-scale energies; de-excitation of nearly-degenerate 
states or annihilation of very light DM might change this 
conclusion to some degree). Modeling of reionization may 
pose a more significant challenge for analyses relying on 
low redshifts (z ~ 10); note, however, that the trans- 
parency of the universe at these redshifts means that in 
realistic scenarios (even decay-like scenarios) the bulk of 
the effect on the CMB comes from earlier times. 



A. Brief review of the Fisher matrix 

The degree to which energy deposition is observable in 
the CMB can be captured by the Fisher matrix for energy 
deposition, denoted F e , which is obtained by contracting 
the transfer matrix T (Equation [5]) with the appropriate 
covariance matrix for the C/s (e.g. [37l - l39| ), 



21+1 



( (cf T ) 2 (cJ E ) 2 



TTr<TE 



{cj e y (c EE y 



r<EEr<TE 
(cf E ) 2 + CJ T C EE 



(Fe] 



E 



dCt 
da. 



dC e 
dai 



(3) 



3 Log binning can of course also be employed for the annihilation- 
like case; there is no clear best choice there, so we will use linear 
binning as the default but show results for both options. See 
Appendix 1X1 for a discussion. 
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For experiments other than the perfect cosmic vari- 
ance limited (CVL) case, noise is included by replacing 

n TT,EE n TT,E. 

noise power spectrum and is given by 



jyTT,EE w Y ieie jy j s t ne e ff ec tive 



(4) 



Here 9 describes the beam width (FWHM = 6>V81n2), 
and the raw sensitivity is (wp)" 1 = (AT x FWHM) 2 , 
with all angles in radians. The standard deviation of 
the parameter ctj, marginalized over uncertainties in the 
other parameters, is given by a ai > (F~ 1 )^ 2 . The pa- 
rameter o,i is then detectable at la if its signal-to-noise 
a.i/<j ai is larger than 1. 

So far, we have not taken into account covariance be- 
tween the standard cosmological parameters and the en- 
ergy deposition parameters, but in fact there are signif- 
icant degeneracies between them. In particular, shift- 
ing the primordial scalar spectral index n, can absorb 
much of the effect of energy deposition [H, HJ. There- 
fore we must marginalize over the cosmological param- 
eters, since the naively most measurable energy deposi- 
tion history may be strongly degenerate with them and 
thus difficult to constrain. We parameterize the usual 
six-dimensional cosmological parameter space by the fol- 
lowing set of parameters: the physical baryon density, 
ujb = flbh 2 , the physical CDM density, ui c = £l c h 2 , the 
primordial scalar spectral index, n s , the normalization, 
A s (k = 0.002/Mpc), the optical depth to reionization, r, 
and the Hubble parameter Hq. 

4 See Appendix[B]for a detailed explanation of this projection and 
how it relates to the standard marginalization prescription. 



Using exactly the same machinery as described above 
for the energy deposition histories, we determine the 
derivatives of the C/s with respect to changes in the cos- 
mological parameters, again assuming that these changes 
are in the linear regime. Then these Ci derivatives are 
vectors spanning an n c -dimensional subspace of the space 
of all Ci derivatives (where for the standard parameter 
set n c = 6); only directions orthogonal to this subspace 
can be constrained. We can regard marginalization over 
the cosmological parameters as simply projecting out the 
components of the energy deposition derivatives orthog- 
onal to this subspace 4 . 

In analogy with Equation[3l we now use the derivatives 
with respect to both energy deposition and the cosmo- 
logical parameters to construct the full Fisher matrix, 



F n = 



F F 

FT F c 



(5) 



where F e is the Fisher matrix for solely the energy de- 
position parameters, F c is the Fisher matrix of the cos- 
mological parameters, and F v contains the cross terms. 
The usual prescription for marginalization is to invert the 
Fisher matrix, remove the rows and columns correspond- 
ing to the cosmological parameters, and invert the result- 
ing submatrix to obtain the marginalized Fisher matrix 
F (e.g. [111]). When the number of energy deposition 
parameters is much greater than the number of cosmo- 
logical parameters, it is convenient to take advantage of 
the block-matrix inversion, 



Fn 



-F~ 1 F'F 



F e — F V F C F v J 



F: 



{F e -F V F-^FJ) 1 F V F- 



1 + F£ (F e — F v F ( T 1 Fj) ^F V F ( 



(6) 



We can now read off the marginalized Fisher matrix as B. Experimental parameters 

F = F e — F V F ( T 1 F^ (note that F has the same units as 

F e ). 



The Fisher matrix approach to estimate detectabil- 
ity is optimistic in the sense that it assumes the likeli- 
hood function is Gaussian about its maximum; for non- 
Gaussian likelihoods, the significance of a given energy 
deposition history will generally be smaller, and any con- 
straints on the amount of energy deposition will be weak- 
ened . We verify that the Fisher matrix method gives 
results consistent with previous studies of WMAP limits 
on constant p a nn in flU CI 



For comparison to the existing literature and con- 
straint forecasting, we consider the WMAP 5, WMAP 7 
and Planck experiments, as well as a theoretical experi- 
ment that is CVL up to I = 2500. The beam width and 
sensitivity parameters for WMAP and Planck are given 
in Table H We use only the W band for WMAP and 
the 143 GHz band for Planck, under the conservative 
assumption that the other bands will be used to remove 
systematics. The effect of partial sky coverage is included 
by dividing E e by / sky = 0.65. 
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Experiment 


Beam 


10 b AT/T W b AT/T 




FWHM (arcmin) 


(I) 


(Q,U) 


WMAP (5 yr, Q band) 


29 


6.7 


9.5 


WMAP (5 yr, V band) 


20 


7.9 


11.1 


WMAP (5 yr, W band) 


13 


7.6 


10.7 


Planck (100 GHz) 


10 


2.5 


4.0 


Planck (143 GHz) 


7.1 


2.2 


4.2 


Planck (217 GHz) 


5.0 


4.8 


9.8 



TABLE I: Detector sensitivities and beams for dif- 
ferent CMB temperature and polarization experiments. 
Results for WMAP temperature sensitivity are t aken 
from 40], with the noise reduced by \/5/4 (\/7/4 
for WMAP 7) to account for the longer integration 
time. The polarization noise for WMAP is taken to be 
v2x the temperature noise. WMAP beam widths are 
taken from [4l]]. The sensitivity and beam width for 
Planck are taken from the Planck Blue Book, available at 
http://www.rssd.esa.int/SA/PLANCK/docs, and assume 14 
months of Planck data. 
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FIG. 2: The effect of the number of included £'s, and the 
number of included frequency bands, on the constraint on 
a constant-p ann energy deposition history; here we show the 
value of Pann corresponding to a 2a signal. 



arbitrary energy deposition history by the "transfer ma- 
trix" T (Equation [5]) we simply compute the effect on 
the Cis of each of a range of non-zero p an n values using 
CosmoRec and CAMB. Since no evidence for energy depo- 
sition has been found to date, we use the predictions of 
the best-fit standard cosmological model (as determined 
by WMAP 7 [H) as a proxy for the WMAP data. Us- 
ing the covariance matrix and prescription for marginal- 
ization laid out in ijll A[ we compute the A% 2 with the 
null hypothesis for each p ann , and by interpolation deter- 
mine how large a p ann would be disfavored at 2a. Our 
results are shown in Figure [5J for both WMAP 5 and 
WMAP 7 noise parameters, with use of one or three fre- 
quency bands (the former is our standard conservative 
approach, the latter is more optimistic), as a function of 
the maximum I included in the analysis (£ m i n is always 
2)- 

In earlier work, [9| found a WMAP 5 limit of p ann < 2 x 
1(T 6 m 3 /s/kg = 3.6 x 1CT 27 cm 3 /s/GeV, using WMAP 5 
data and the CosmoMC code to perform a full likelihood 
analysis, in reasonable agreement with our estimate for 
the WMAP 5 1-band case. Using WMAP 7 data, 0] 
provided an updated constraint p ann < 2.4 x 10~ 27 
cm 3 /s/GeV, also in very good agreement with our es- 
timate. Note that [t| [lj] included Lyman-a photons in 
their analysis, which we have not done here and which 
would strengthen the constraints slightly: the small dif- 
ferences between our constraints and those in the litera- 
ture can probably be ascribed to the combination of this 
factor, our use of CosmoRec rather than RECFAST, the fact 
that we use the best-fit standard cosmological model as a 
proxy for the real WMAP data, and the optimism inher- 
ent in a Fisher matrix analysis. It is reassuring that none 
of these factors seem to give rise to large discrepancies in 
the results. 

We see that the WMAP limits are essentially un- 
affected by £ max for 4nax <; 1000, and the projected 
Planck bound appears stable for £ maK > 1500. For the 
CVL case, of course, higher Ps will always yield more 
information, but the rate of improvement with £ max is 
quite slow for the £'s we are considering. 



D. Numerical stability of derivatives and linearity 



C. Comparison to previous results 

Constraints on energy deposition from WMAP 7 have 
been studied previously in the case where p ann (z) = p aim 
is constant [Uj]. We have obtained analogous con- 
straints using the estimates for experimental sensitivity 
described in mi Atill Bl taking the six standard cosmo- 
logical parameters and p ann as our parameter set, and 
marginalizing over the cosmological parameters. 

Note that in this simplified scenario there is no need 
for the Gaussian basis functions described at the start 
of the section, and rather than describe the effect of an 



When dealing with general energy deposition histories, 
we hope to work in a regime where the effect of deposi- 
tion on the CMB is linear, so that the effect of a gen- 
eral energy deposition history can be described in terms 
of a linear combination of basis energy deposition his- 
tories. This is the idea behind characterizing the effect 
of new parameters entirely in terms of the transfer ma- 
trix of derivatives, T, and the Fisher matrix F derived 
from it. Equivalently, linearity means it is sensible to 
speak of a single transfer matrix T largely independent 
of the "fiducial" energy deposition history about which 
the derivatives dCifdon are taken (our default assump- 
tion is that this "fiducial" energy deposition is zero). If 
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the energy deposition history being studied is too great a 
perturbation away from the fiducial, the first derivatives 
will no longer accurately describe its effect on the Ce's, 
and the Fisher matrix estimate of its significance will 
break down. In this subsection we discuss the numerical 
stability of the derivatives, and the degree to which they 
describe the effect of arbitrary energy deposition histories 
on the Ci's. 

To obtain derivatives of the Ci spectra, we calculate a 
grid of Ci values for each of the parameters considered 
in the Fisher matrix. For each i and each parameter, the 
derivative is extracted from a polynomial fit to Ci as a 
function of the parameter values. 

We estimate that we have calculated derivatives to 
a precision of ~ 1% for cosmological parameters. For 
energy deposition, the error is less than ~ 2% for the 
most relevant redshifts of z < 800 and rises to ~ 5% 
for z < 1100. While the error in the derivatives becomes 
larger for higher redshifts, the effect on the PCA is small; 
the difference in signal-to-noise is at the percent level for 
WMAP 7 and at the few percent level for Planck and 
the CVL case. The dominant numerical error here comes 
from the numerical accuracy limitations of CosmoRec and 
CAMB. We estimate the error by comparing the derivatives 
obtained from two different accuracy settings 5 . 

The derivatives used in the Fisher matrix are evaluated 
at the fiducial cosmology (with no energy deposition). 
The assumption of linearity is that these derivatives are 
still correct away from the fiducial. For the standard set 
of six cosmological parameters, the biases to the cosmo- 
logical parameters induced by the maximum permitted 
energy deposition from WMAP 5 generally lie well within 
the linear regime. 

For large energy deposition, the effect on the CV's 
is nonlinear, i.e. not directly proportional to the de- 
posited power as parameterized by the a*; equivalently, 
the derivatives about a fiducial large energy deposition 
are not the same as for zero energy deposition. Our poly- 
nomial fits for the derivatives, described above, also allow 
us to check the extent to which nonlinearity may become 
important: that is, the extent to which O{o%) corrections 
to the effect on the Ce's are non- negligible. 

The amount of energy deposition such that nonlinear- 
ities become important depends on redshift z. This can 
be estimated by the fractional rate of ionization per Hub- 
ble time, (dni on /dt)/(n® on H(z)), arising from the energy 
deposition (where dn ion /dt is related to dE/dtdV accord- 
ing to the prescription of [42]). For two fiducial cases this 
quantity is shown in Figure [T] Conversely, the energy de- 
position at redshift z such that {drn on / dt) / (n® on H (z)) 
1 gives a measure of what energy deposition is required 



1.000 



5 The two settings we use are (1) runmode=0 and cosmorecacc=0 
and (2) runmode=l and cosmorecacc=2. In general, we use CAMB 
accuracy settings accuracy_boost=4, l_accuracy_boost=4, 
l_sample_boost=4 for computing derivatives. For the MCMC 
runs, we set all three parameters to 1 for speed. 
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FIG. 3: The degree of nonlinearity in the computed signifi- 
cance of a sample energy deposition history, for p an n constant, 
using WMAP 7 noise parameters. We show the ratio of (1) 
the S/N estimated by a linear extrapolation from small en- 
ergy deposition to (2) the "true" S/N (estimated as in C|l . 
as a function of p ann . The solid, dashed and dotted lines in- 
dicate the WMAP 7 2a upper limit on p ann , the value of p ann 
for which the nonlinearity is 10%, and the value for which the 
nonlinearity is 1%, respectively. The red dot-dashed line in- 
dicates the 2a upper limit on p ann that would be obtained by 
linearly extrapolating the significance from small energy de- 
position, which overestimates the significance and hence leads 
to a too-strong constraint. 



before nonlinearities may become significant. For each 
redshift bin, we use the polynomial fits of 5Ce{on) to nu- 
merically calculate the derivatives at this level of energy 
deposition. We then find 1% corrections (averaged over 
i) to the fiducial derivative (pCej 'dati)\ ai =o. 

One simple test of the effect of nonlinearity is the de- 
gree to which the "true" bound on constant p ann from a 
given experiment, evaluated without assuming linearity 
as described in $11 C[ differs from the bound we would 
obtain by taking derivatives dCej dp ami at p ann = 0, and 
assuming linearity, i.e. taking the effect on the Ce's to be 
given by (dCe/dp ann ) | Pann =o- Equivalently, we can 
compare the signal-to-noise estimate for the two meth- 
ods, as a function of p ann . 

In Figure [3] we show an example of this test using 
WMAP 7 noise parameters. For each value of p ann , we 
compute both the S/N of the resulting signal directly (as 
in mi C|) , and the significance that would be obtained by 
a linear extrapolation from small p ann - The ratio of the 
extrapolated S/N to the true S/N is 1 + e, and e provides 
a measure of nonlinearity. We see that e approaches 0.3 
close to the 2a upper bound from WMAP 7, so the esti- 
mated 95% bound from WMAP 7 is roughly 30% weaker 
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than would be expected from a linear extrapolation from 
small energy deposition, but the degree of nonlinearity 
falls rapidly at lower energy deposition. We will see later 
that for the energy deposition probed by Planck, e is 
only a few percent. Note that even for WMAP 7, the 
difference in the bound is almost entirely due to an over- 
all ^-independent normalization factor; the shape of the 
SCe's is almost unchanged. 

III. PRINCIPAL COMPONENT ANALYSIS 

The effects of energy deposition at different redshifts 
on the Cgs are highly correlated, and so the effects of a 
large class of energy deposition histories can be charac- 
terized by a small number of parameters. Principal com- 
ponent analysis provides a convenient basis into which 
energy deposition histories can be decomposed, with the 
later terms in the decomposition contributing almost 
nothing to the effect on the Cis. It thus allows gen- 
eralization of constraints on energy deposition to a wide 
range of models (subject to the linearity assumption dis- 
cussed above). 

A. The principal components 

Having obtained the marginalized Fisher matrix F, di- 
agonalizing F: 

F = W T AW, A = diag(Ai,A 2 ,....,A N ) (7) 

yields a convenient basis of eigenvectors or "principal 
components" . W is an orthogonal matrix in which the 
i-th row contains the eigenvector corresponding to the 
eigenvalue A^. If we compute derivatives for N redshift 
bins, then the N x N Fisher matrix has N principal com- 
ponents. The eigenvectors are orthonormal in the space 
of vectors i — 1..N. Let us label these vectors e», 

with corresponding eigenvalues Aj, i = 1..N. Our con- 
vention is to rank the principal components by decreasing 
eigenvalue, such that e\ has the largest eigenvalue. 

Note that the principal components may be signifi- 
cantly different from the unmarginalized principal com- 
ponents, or the eigenvectors of F e . Figure @] shows 
the first three principal components for WMAP 7, 
Planck and a CVL experiment, both before and after 
marginalization, for the annihilation-like case (dE/dt oc 
Pann(z)(l + z) 6 ) with 50 linearly-spaced redshift bins. We 
see that while the shapes of the PCs are qualitatively sim- 
ilar, marginalization produces noticeable changes to the 
PCs, as does changing from one set of experiment pa- 
rameters to another. The differences become more pro- 
nounced for higher PCs. 

Note that the shapes of the principal components can 
be affected by a number of other different factors: choice 
of binning, choice of ionization history calculator, en- 
ergy deposition model, fiducial cosmological model con- 
sidered, etc. We discuss these effects in Appendix |A"1 



In Figure [5] we show the first six marginalized PCs for 
Planck, for annihilation-like (dE/dt oc p an n(z)(l + z) 6 ) 
and decay-like (dE/dt oc p,\ ec (z)(l + z) 3 ) energy deposi- 
tion histories. We show the annihilation-like case with 
both log and linear binning. We note that the first 
principal component is always largely or completely non- 
negative, and (in the annihilating case) peaked around 
redshift 600. The first PC can be thought of as a weight- 
ing function, describing the effect of energy deposition on 
the CMB (orthogonal to the effect of shifting the cosmo- 
logical parameters), as a function of redshift 6 . 

In Figure [5] we show the effect on the ionization his- 
tory for the first three Planck PCs in the annihilation 
case, with each PC multiplied by an energy deposition 
coefficient of e = 2 x I0~ 27 cm 3 /s/GeV to obtain p ann (z). 
Note that this energy deposition is too large to be strictly 
in the linear regime; this figure illustrates the shape and 
size of the effect in the linear regime, the true effect for 
this value of e will be somewhat smaller. 

As previously, we have considered "annihilation-like" 
and "decay-like" energy deposition histories separately. 
If both analyses were performed over the same redshift 
range, then while the principal components might ap- 
pear different, they would span the same space of energy 
deposition histories. If all principal components were re- 
tained, the difference between the two would simply be 
equivalent to a change of basis, and provided sufficient 
principal components are retained, this will still be ap- 
proximately true. However, a particular energy depo- 
sition history may be described by the early principal 
components much better in one case than in the other; 
in particular, energy deposition histories for which the 
effect on the CMB is dominated by low redshifts will 
not be well described by the (first few of the) default 
annihilation-like PCs. Thus we present results for both 
cases. 

B. Mapping into SC'e space 

Let us consider the mapping into 6Ci space of these 
marginalized principal components. Applying the trans- 
fer matrix T (Equation [5]) to the eigenvectors yields 
a set of N vectors in the space of Ce perturbations, 
5Ce = Tei = hi. The h^s should be understood as SC^s 
per energy deposition, and have units of Ci/p(z). 

We can define a dot product on the space of SCe's by 

h * ■ h o ■ = /l « E 7 1/l ^ = e T F * e 3 ( 8 ) 
t 

We then see that while the PCs are orthogonal, the hi are 
in general not orthogonal to each other, nor to the SCe's 



Note that the shift in the peak position between log and linear 
binning is to be expected, as one "weighting function" would be 
integrated over dz and the other over ciln(l + z)\ see Appendix 
|X]for further discussion. 
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FIG. 4: The first three principal components for WMAP 7, Planck and a CVL experiment, both before and after marginalization 
over the cosmological parameters. 




FIG. 5: The first six principal components for Planck after marginalization, in the case of (left) annihilation-like redshift 
dependence with linear binning, (center) annihilation- like redshift dependence with log binning, and (right) decay- like redshift 
dependence with log binning. Note that for decay-like energy deposition histories, the redshift range is extended down to z = 10 
in order to fully capture the effect on the CMB - see [jll] This larger redshift range makes linear binning impractical. 



from the cosmological parameters. They correspond to 
actual energy deposition histories, and in general, there 
is no such history that is precisely orthogonal to all the 
cosmological parameters. 

However, we may decompose the hi into components 
parallel and perpendicular to the space spanned by vary- 
ing the cosmological parameters, and denote the perpen- 
dicular components h^. The projection operator that im- 
plements this decomposition is given in Appendix[Bj The 
vectors are now orthogonal amongst themselves, as 
well as to the cosmological parameters, and their norms 
are given by the square root of the marginalized eigenval- 
ues Xi. It is these hj~'s which determine the detectabil- 
ity of the marginalized principal components, and which 
form an orthogonal basis for residuals which cannot be 
absorbed by varying the cosmological parameters. The 
addition of the parallel components, to recover the /ij's 
from the hj-'s, ensures that the /i^'s correspond to energy 
deposition histories, and so provide an orthogonal basis 



in redshift space. 

In Figure [3 we show the mapping of the first three 
(marginalized) PCs for Planck into the space of r5CVs; in 
Figure [51 we show the components of these 8CiS which 
are orthogonal to the space spanned by varying the cos- 
mological parameters. Figure |H] demonstrates this pro- 
jection for a sample DM annihilation model, summing 
over principle components, and decomposing the effect 
on the C/s into components perpendicular and parallel 
to the cosmological parameters. 

The eigenvectors of the Fisher matrix {e^} thus provide 
an orthogonal basis in both relevant spaces, and their 
eigenvalues precisely describe the measurability of a "unit 
norm" energy deposition history with z-dependence given 
by the eigenvector. For an arbitrary energy deposition 
history which we now write as 

N 

p(z) =^2e l e l (z), (9) 

z=l 
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error bars, and will be challenging to measure or con- 
strain. 

We now outline the method for reconstructing and con- 
straining the PC coefficients, or any specific energy depo- 
sition history, using the PCA formalism. We investigate 
the number of PCs that can generically be measured at 
> 1(7 by Planck and a CVL experiment, for arbitrary 
energy deposition, and show results for broad classes of 
example models. We also consider the biases to the cos- 
mological parameters that are induced if energy deposi- 
tion is present but ignored; we present results for each 
principal component, so the biases due to an arbitrary 
energy deposition history can be immediately calculated. 
Our estimates of detectability and the biases will be ver- 
ified using CosmoMC in the 9 VII 



FIG. 6: Fractional change to the ionization fraction x e in the 
presence of energy deposition, for the first three (marginal- 
ized) principal components in Planck. The curve shown is ex- 
trapolated from the linear (small energy deposition) regime, 
with normalization factor £1,2,3 = 2 x 10 -27 cm 3 /s/GeV. 



the expected A^ 2 relative to the null hypothesis of no en- 
ergy deposition is £ 2 Ai. If the Si coefficients are com- 
parable, the relative sizes of the eigenvalues describe the 
fractional variance attributable to each principal compo- 
nent (eigenvector). 

A brief comment on unit conventions: we take the 
{ei} and {Gi} to be dimensionless, with the units of 
p(z) (cm 3 /s/GeV) carried by the coefficients cti, Ei. The 
derivatives (and transfer matrix) then have units of 
Ci/p(z), and the Fisher matrix and its eigenvalues have 
units of l/p(z) 2 (since the covariance matrix E has units 
of C\ ) . Note also that due to the units of the covariance 
matrix, the dot product defined above takes two vectors 
in Ci -space to a dimensionless number (if the vectors 
have units of Ci). 



IV. DETECTABILITY 

For a general energy deposition history, the PCs pro- 
vide a basis in which, by construction, the basis vectors 
are ranked by the significance of their effect on the C/s. 
The measurability of a generic (smooth, non-negative) 
energy deposition history can thus be accurately de- 
scribed by the first few PCs 7 . Equivalently, the coeffi- 
cients of later principal components have extremely large 



7 It is in principle possible for the coefficients Ei to be zero for 
i < n for some n, but if n is large this implies a very unphysical 
energy deposition history that oscillates rapidly between positive 
and negative values. While "negative energy deposition" might 
perhaps have a physical interpretation in terms of increased ab- 
sorption of free electrons, such an interpretation is not at all 
obvious, and so we focus on smooth, non-negative energy depo- 
sition histories. 



A. Estimating limits from the Fisher matrix 

As mentioned previously, the perpendicular compo- 
nents of the SCg's, hj~, are orthogonal with norms \fXi. 
They are also orthogonal to the space spanned by vary- 
ing the cosmological parameters. Given these results 
and a measurement of the temperature and polarization 
anisotropics, it is straightforward to estimate general con- 
straints on the energy deposition history from the Fisher 
matrix formalism. Note that in a careful study, one would 
instead use CosmoMC to perform a full likelihood analysis, 
using the Fisher matrix results only to determine the op- 
timal principal components, as we demonstrate in flVIl 
We outline this simple method only to help build intu- 
ition and to clarify later comparisons between the Fisher 
matrix method and the CosmoMC results. 

The first step is to extract any residual between the 
data and the best-fit model using the standard cosmologi- 
cal parameters; let us denote this residual by Rj T ' EE TE . 
Then we take the dot product (as defined in Equation [8]) 
of this residual with the hf~ vectors, normalizing by the 
corresponding eigenvalues (this normalization is required 
because the h^s are orthogonal, but not orthonormal; 
see Appendix [Bjl : 



Xi 



(10) 



The resulting are the model-independent reconstructed 
coefficients for the marginalized principal components. In 
the absence of energy deposition, we expect them to be 
zero (within uncertainties). 

The individual la uncertainties on each of these co- 
efficients are 1/y/Xl, in the sense that if a single coeffi- 
cient is perturbed away from its best- fit value by 1/ y/Xl, 
the corresponding energy deposition history will be dis- 
favored at la. Thus it is possible to set a very general 
model-independent constraint on each of the coefficients, 
Ei = e l ± -j= (at la). 

Given an arbitrary energy deposition history, we can 
decompose it into the principal components, each with its 
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FIG. 7: The mapping of the first three principal components for Planck, after marginalization, into SCe space. The PCs are 
multiplied by £i(z) — 2 x 10~ 27 cm 3 /s/GeV for all i, to fix the normalization of the <5CV's. 




FIG. 8: The _L components of the first three principal components for Planck, after marginalization, mapped into SCe space. 
The normalization is the same as for Figure [7] 



own coefficient, and compare these coefficients £j to the 
bounds. For any particular model, a stronger constraint 
can be set by noting that, 

A x 2 =Y / K(e l -e l f. (11) 

i 

This A% 2 is relative to the best-fit model including both 
energy deposition and the standard cosmological param- 
eters; the A\ 2 relative to the best-fit standard cosmolog- 
ical model 8 is simply J2 i Ai£j(£, — 2Si). 

This method has the usual deficiencies of the Fisher 
matrix approach: it assumes a Gaussian likelihood and 
also linearity of the derivatives, and so can only be used 
for an estimate. In WII we will go beyond the Fisher 
matrix approach and present constraints derived from a 
likelihood analysis using CosmoMC: in the same way as 
this estimate, those limits can be expressed as bounds on 



Of course, if the best-fit energy deposition history is everywhere 
zero, i.e e; Rj for all i, these two quantities are identical. 



(a simple combination of) the PC coefficients, and will 
therefore be immediately applicable to a wide range of 
models for energy deposition. 

B. Sensitivity of future experiments 

For an energy deposition history where the sizes of the 
coefficients, | j , are all similar, the respective detectabil- 
ity of the PCs are given simply by their eigenvalues. Lit- 
erally taking all the coefficients to be the same docs not 
give a physical energy deposition history (since the later 
eigenvectors are highly oscillatory) , but it is in some sense 
a "generic" scenario: none of the PCs have coefficients 
that are fine-tuned to be small, so slight changes to p(z) 
or the basis of PCs are unlikely to drastically change the 
detectability of the different components. 

We define detectability of the PCs with respect to this 
"generic" case; of course, detectability of any particular 
model depends on the relative sizes of coefficients. We 
consider a number of physical examples below to illus- 
trate that, in some sense, the generic case is a reasonable 
average over a wide class of models of interest. 
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FIG. 9: Decomposition of SCe from energy deposition with constant p a nn(z) into parallel (||) components which can be absorbed 
by changes in the cosmological parameters, and perpendicular (_L) components that cannot be absorbed by such changes. The 
overall effect of the energy deposition is suppression of high-£ modes, due to the increased optical depth, and enhancement of 
low-^ polarization modes, as discussed in The suppression at high £ is clearly seen in the TT and EE spectra; the effect is 
also present in the TE spectra, with the peaks of SCf E occurring at the troughs of Cj , and vice versa. The normalization 
here is p ann = 2 x 10 -27 cm 3 /s/GeV, comparable to the latest limits from WMAP 7+ACT [lj. This decomposition depends 
on the sensitivity of the experiment; the case shown is WMAP 7 single band. 



As discussed previously, [10| derived a set of energy 
deposition profiles corresponding to a range of DM an- 
nihilation models. These models provide a convenient 
set of example energy deposition histories, although they 
all have very similar effects on the CMB (see SjV|). We 
adapt the code developed in [l(| and discussed in detail 
there to obtain similar physical f(z) curves for the case 
of decaying dark matter with a long lifetime. 

While the DM itself must have a lifetime consider- 
ably longer than the age of the universe, there could be 
other metastable species which decay during the redshift 
range we study (z ~ 10 — 1300), or excited states of the 
dark matter which decay to the ground state + Standard 
Model particles (e.g. |4^I447| and references therein). In 
this case the decay rate would cut off exponentially for 
z < z(t), although heating and ionization of the gas 
could continue for some time after that: we can again 
obtain detailed p(z) curves for different decay lifetimes 
using the methods of Models of this type provide a 
simple class of examples suitable for use with the PCs de- 
rived for the case of decaying DM, since the underlying 
dE/dt cx (1 + z) 3 redshift dependence is the same (al- 
though for models with lifetimes short enough that the 
energy deposition has ceased shortly after recombination, 
the PCs derived for the annihilation-like case may work 
better). 

For the annihilating case, asymmetric dark matter sce- 
narios can furnish a similar set of examples [481 - 4531 . In 
such scenarios the DM sector possesses an asymmetry 
analogous to that in the baryon sector, and it is this 
asymmetry which sets the DM relic density rather than 
the annihilation cross section. In the minimal case there 
is thus no requirement for an annihilation signal in the 
present day or during the epoch of recombination, but it 
is nonetheless possible to have a large late-time annihi- 



lation signal, by repopulation of the depleted component 
at late times, or by oscillations from the more-abundant 
to the less-abundant component [53|, [HJ . As a simple 
example, we consider models where another species de- 
cays to repopulate the less-abundant DM state [53j , thus 
causing the annihilation to "switch on" as 1 — e _t / r at a 
characteristic timescale r (with z T being the correspond- 
ing redshift). We compute the p(z) curves for a range of 
r. Finally, for both annihilation and decay we consider 
the constant p(z) case, studied in [36[ (for decay) and 
H Q (for annihilation) , to facilitate comparison with the 
literature. 

Figure [TU] shows the detectability of the principal com- 
ponents in Planck and the ideal CVL experiment for 
these annihilating and decaying models, with the en- 
ergy deposition normalized to lie at the 95% limit from 
WMAP 7. In the "generic" case, we set the sizes of the 
coefficients of the Planck (or CVL) PCs to be \ei\ = e = 

MAP - The actual WMAP 7 signal-to-noise 
for the model is 



S_ 

N 



WMAP 



Planck WMAP 



1/2 



and thus depends on the signs of £j, but the generic case 
is meant to indicate the typical detectability for a class 
of models, so we instead use the WMAP 7 constraints to 
set an overall scale for \si\. 

We also show the detectability for each PC if p(z) cx 
&i{z), or assuming the energy deposition history has zero 
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FIG. 10: (a) The sensitivity for Planck (single-band), after marginalization, for various models subject to constraints from 
WMAP 7 single-band at 2a. The left figure assumes annihilation-like energy deposition and the right figure assumes decay-like 
energy deposition. The top panels show: (1) assuming p(z) oc ei{z) for each PC, (2) the generic case where all PC coefficients 

have equal magnitudes [sj| = e — 2/y^5]V \ f^ MAP , (3) constant p(z), and (4) taking p(z) oc f(z), with f(z) from the models 

in [l(J. For the left figure, the hatched region indicates the range of results from changing the ionization history calculator 
and including or neglecting the effects of helium and Lyman-a photons, described in £ ]A 31 A 4l The bottom panels show some 
sample z T models for asymmetric annihilating dark matter (left) and decaying species (right), as discussed in ^IV Bl (the labels 
describe the initial particle mass, and the SM final state for annihilation or decay), and an extreme case where p(z) = for 
200 < z < 900 and constant outside that range, (b) Same as (a), but for a CVL experiment, (c) The models in (a) and (b) 
for annihilation-like (left) and decay-like (right). 
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overlap with all other PCs 9 . As mentioned previously, 
this is not a physical assumption (requiring an "energy 
deposition" oscillating rapidly between positive and neg- 
ative values): in such a case the effect on the C/s is so 
small that the normalization of the "energy deposition" 
could be very large and still consistent with WMAP. Con- 
sequently, arbitrarily high PCs can be measured if they 
are the sole contributors to the energy deposition history. 

We see that models with decay-like redshift depen- 
dence and those with annihilation-like redshift depen- 
dence tend to have roughly the same number of measur- 
able parameters. In both cases, generally 2-3 components 
are potentially measurable in Planck and up to 5-7 for a 
CVL experiment. 

As a side note, the improvement of these constraints 
between WMAP 7 and future experiments is in large part 
due to (anticipated) better measurements of the polariza- 
tion. In the absence of polarization data (i.e. using the 
TT spectrum only), we would expect the constraints to 
weaken by a factor of ~ 3 for WMAP 7,^7 for Planck , 
and ~ 14 for a CVL experiment. Here we have taken 
the square root of the eigenvalue of the first principal 
component as a proxy for sensitivity, which will be ap- 
proximately true for models with a non-negligible overlap 
with the first PC. 



C. Biases to the cosmological parameters 

If energy deposition is present but neglected, it can 
bias the measurement of the cosmological parameters by 
a significant amount. For WMAP, the partial degeneracy 
between varying n s and the effects of energy deposition 
means that the dominant bias is a lcr negative shift to n s . 
The improved polarization sensitivity of Planck largely 
lifts the degeneracy with n s , but due to the smaller error 
bars of Planck other parameters develop non-negligible 
biases: at the maximum energy deposition allowed by 
WMAP 7 at 2cr, Planck parameter estimates are generi- 
cally biased at > la for w c , H , and A s . 

Calculation of the biases is exactly complementary to 
calculating the marginalized Fisher matrix. While the 
marginalization can be understood as projecting out the 
degeneracies with the cosmological parameters, the bi- 



9 If the PCs were the same for the different experiments, this would 
give an upper bound on the detectability of the ith PC, given 
WMAP 7 2<r constraints. However, the PCs for different exper- 
iments are not orthogonal, e f lanck ■ e^MAP str j ct 
upper bound for the S/N of the ith Planck PC is given by 

(S/N) Planck < 2\J \ planck y \ £ Planck . £ WMAP j I yW M AP | 

with the analogous result for a CVL experiment. However, this 
quantity is not very useful as an upper bound; for example, if p(z) 
is proportional to a high WMAP PC, the normalization of p(z) 
is essentially unconstrained, but the detectability for Planck may 
be very significant if there is even a small overlap with the first 
Planck PC. 



ases are given precisely by the effect of energy deposition 
in those degenerate directions. To be precise, suppose 
that some eigenvector ej has true coefficient ej ^ and 
we falsely assume Ej to be zero: then each of the cosmo- 
logical parameters 9i will be shifted by an amount S9i. 
The matrix of derivatives dO.Jdej, i = l..n c , j = 1..N, is 
given simply by J2k { F c lF I) ik (ej)k- 

Thus we can partition the biases into the bias per PC, 
which is shown in Figure [TT] for WMAP 7 and Planck. 
For a generic energy deposition history, the total bias is 
dominated by the bias from the first few PCs, consistent 
with the fact that later PCs are undetectable and can es- 
sentially be neglected in any fit to the data. As expected 
from [12|, the largest bias for WMAP 7 is to n s . 



V. A UNIVERSAL Pa . nn (z) FOR WIMP 
ANNIHILATION 

Solutions for the redshift dependence of the efficiency 
function f(z) (and hence the energy deposition history 
Pann(z)), for 41 different combinations of dark matter 
mass and annihilation channel, were presented in [loj . 
We can use these 41 energy deposition histories, rather 
than 5-functions in z, as the input states for a principal 
component analysis, specialized for the particular case of 
conventional WIMP annihilation. Even after marginal- 
ization over the other cosmological parameters, we find 
that in this case the first eigenvalue completely dominates 
the later ones, accounting for 99.97% of the total vari- 
ance in WMAP 7, Planck and the CVL forecast: thus, 
to a very good approximation, for any of the DM mod- 
els studied in [10| (or any linear combination of the final 
states studied there) , the effect on the Cg 's is determined 
entirely by the dot product of p ann (z) with the first PC, 
with the ^-dependence given by mapping the first PC to 
CVspace. 

This conclusion agrees with the statements in [l^, EH 
that the effect of DM annihilation can be captured by a 
single parameter. To put it another way, given equal co- 
efficients for the first two principal components (which is 
already rather conservative, since the f(z) curves studied 
are generally very similar to the first PC), the signal cor- 
responding to the first PC would be roughly 60 x larger 
than the signal corresponding to the second PC: the ex- 
istence of energy deposition would have to be detected at 
60(7 for even a lcr measurement of the second component 
to be possible. Measurements of the later components 
would be far more difficult still: the sum of the first two 
eigenvalues accounts for 1 — 4.8 x 10~ 7 , 1 — 6.0 x 10~ 7 , 
and 1 - 8.3 x 10" 7 the total variance in the WMAP 7. 
Planck and CVL cases respectively. The effective /-value 
of various WIMP annihilation models is then just given 
by the dot product of their f(z) curves with this first 
principal component. We provide effective /-values for 
all models considered in [lOj on our website. 

Let us denote this first principal component by 
ewiMp(z). We note that the differences between the 
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FIG. 11: For the iih PC, the contribution to the bias to cosmological parameters in WMAP 7 (left panel) and Planck (right 
panel), relative to the error bars forecast from the Fisher matrix. The normalization is that of the "generic" case (see discussion 
in i]IVBI or Figure [TO)) . where each PC coefficient has the same absolute value and the overall normalization is the maximum 
allowed by WMAP 7 at 2a. The total bias for the parameter 6 is £L S9i. 



ewiMp(z) curves corresponding to different experiments 
(WMAP , Planck and a CVL experiment) are extremely 
small, at the sub-percent level for all redshifts, and so it is 
reasonable to speak of a single such curve. We can choose 
to normalize ewiMp(-z) so that when it is multiplied by a 
given (dimensionful) factor e to obtain an energy deposi- 
tion history p ann (z) = eewiMp(z), the significance of the 
resulting signal in a particular experiment is the same as 
that of a constant-paim energy deposition with p ann = e. 
This normalization simplifies comparisons with the ear- 
lier literature, in that the constraint on e is precisely the 
same as the familiar limit on constant p a nn- The choice 
of experiment affects the normalization at the percent 
level; as a default, we will use the normalization appro- 
priate for Planck. We plot the resulting e\yiMp(z) curve 
in Figure [12] this curve and the corresponding Ci shifts 
are also available online in tabulated form (see Appendix 

We could of course also write ewiMp(z) m terms of 
the previously calculated principal components, just like 
any other energy deposition history. For example, for 
the Planck PCs derived above, an energy deposition e x 
ewiMp(z) = J2i £ i e i corresponds to {ei, e 2 , £3, ••■} = e x 
{4.64,-0.396,3.11,...}. In particular, this implies that 
the shapes of the <5CVs induced by WIMP annihilation 
are quite close to those shown for the first PC in Figures 

In Figure Q2] we show, for reference, the modifica- 
tion to the ionization history associated with the ewiMP 
curve, and the contributions from the first three principal 
components. We see that the first three PCs provide a 
good description of the ionization history modifications 
at 300 < z < 1000, and in both cases the ionization frac- 
tion is essentially unaffected for z > 1000; we can infer 
that the effect on the CMB of the discrepancy at lower 
redshifts is small. 



VI. COSMOMC RESULTS 

Everything we have done so far assumes both linearity 
and that the Fisher matrix is an adequate description 
of the likelihood function. We now present results of a 
full likelihood analysis using the CosmoMC Markov chain 
Monte Carlo code, in particular examining the biases to 
the cosmological parameters and the detectability of the 
PCs. Throughout this section we use RECFAST 1.5 as our 
ionization history calculator, since we have established 
that the PCs are unaffected by this choice (see £|A 3[) and 
the interface to CosmoMC is better established. 

We sample the six cosmological parameters wj,, ui c , n s , 
lnlO^A^k = 0.002/Mpc), r and H , all with flat pri- 
ors. We consider purely adiabatic initial conditions. We 
parameterize the energy deposition due to dark matter 
annihilation using the marginalized principal components 
in redshift space e, presented in mill We thus include 0, 
1, 3, 5 or 7 additional parameters corresponding to the co- 
efficients of the principal components with highest signif- 
icance, as determined from the Fisher Matrix analysis 10 . 
We impose flat priors on these parameters. Our treat- 
ment of the energy deposition is the same as described 
in the previous sections; we do not include Lyman-a as 
a default (see 4|) . 

The MCMC convergence diagnostic tests are per- 
formed on 4 chains using the Gelman and Rubin (vari- 
ance of chain mean) / (mean of chain variances) R - 1 



We use odd numbers of principal components to illustrate the 
effects of including more PCs because for near-constant p ann , 
the even-numbered principal components tend to have very small 
coefficients, and so including them does not substantially change 
the results. 
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FIG. 12: The universal ewiMP curve (solid black line), nor- 
malized as discussed in the text. Note that the principal 
component analysis is only performed in the redshift range 
between the vertical red solid lines; outside these lines, we 
still plot the linear combination of the input energy deposi- 
tion histories that has been identified as the first principal 
component, to serve as a canonical energy deposition history 
for WIMP annihilation, but the Ct's are not sensitive to the 
details of this energy deposition. In the PCA region we also 
plot the second (dashed) and third (dotted) principal compo- 
nents, with arbitrary normalization. 
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FIG. 13: Fractional change to the ionization fraction x e in 
the presence of energy deposition, for the universal ewiMP 
curve (red, solid). We also write the ewiMp(z) curve as a 
linear combination of the Planck principal components, and 
show (in black) the effect on the ionization history of the 
first three principal components individually (weighted by 
their contribution to ewiMp), and their (weighted) sum. The 
curves shown are extrapolated from the linear (small energy 
deposition) regime, with normalization factor e = 2 x 10 -27 
cm 3 /s/GeV for the ewiMP curve. 
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PC1 x10- 6 PC2 x10- 6 PC3 x10- 6 

FIG. 14: Constraints from the seven-year WMAP data (red), 
and from simulated data for Planck (blue) and a cosmic vari- 
ance limited experiment (green) . The plot shows marginalized 
one-dimensional distributions and two-dimensional 68% and 
95% limits. The mock data for Planck and the CVL experi- 
ment assumed no dark matter annihilation. Three Principal 
Components were used in each run to model the energy de- 
position from dark matter annihilation. The units of the PC 
coefficients here are in m /s/kg, with 1 x 10" 6 m 3 /s/kg = 
1.8 x 1CT 27 cm 3 /s/GeV. 



statistic for each parameter. Our constraints and the 
1 — D and 2 — D likelihood contour plots are obtained 
after marginalization over the remaining nuisance param- 
eters, again using the programs included in the CosmoMC 
package. We use a cosmic age top-hat prior of 10 Gyr 
< t < 20 Gyr. 

We first determine the constraints on parameters using 
the seven-year WMAP data [l| (temperature and polar- 
ization) with the routine for computing the likelihood 
supplied by the WMAP team. For this case only, we 
also marginalize over a possible contamination from a 
Sunyaev-Zeldovich component (see e.g. [Hj]). 

We then generate simulated data for Planck and a CVL 
experiment using a fiducial cosmological model given by 
the best fit WMAP7 model. We simulate the data as- 
suming in one case no energy deposition, and in an- 
other case an energy deposition history with constant 
Pann = 1-78 x 10~ 27 cm 3 /s/GeV. We model experimental 
noise as described in Equation H) For the Planck ex- 
periment, we use the specifications reported in Table U 
for the 147 GHz channel only. For the CVL experiment, 
we limit the maximum resolution of the experiment to 
(max = 2500. In each case, we use the PCs developed for 
that particular experiment, as described in mill 
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Number of 


PC 




WMAP7 


PCs used 






95%c.l. 


1 


PC 


< 1.2 x 10" 


26 cm 3 /s/GeV 


1 


ewiMp(z) 


< 2.43 x 10" 


27 cm 3 /s/GeV 



TABLE II: Upper limits on the first principal component amplitude using WMAP 7 data. In the last line we also show the 
constraints on the ewiMp(z) principal component presented in ijV] The uncertainties reported are upper limits at the 95% c.l. 
We show the constraints obtained when only the first principal component is varied, together with the cosmological parameters. 
We assume in this case a flat positive prior on the amplitude of the principal component. 



Number of 
PCs used 


PC 


Pann = 




Planck 

Pann = 1.78 X 10" 27 cm 3 /s/GeV 


Pann = 




CVL 

Pann = 1.78 x 10~ 27 cm 3 /s/GeV 


1 


PC 1 


(0.2 ± 1.1) X 10" 


-27 


(8.9 ± 1.6) x 10~ 27 


(0.1 ±4.9) x 10" 


-28 


(9.04 ±0.81) x 10~ 27 


3 


PC 1 


(0.4 ± 1.1) x 10" 


-27 


(8.8 ± 1.5) x 10" 27 


(0.6 ±5.0) x 10" 


-28 


(8.71 ±0.81) x 10" 27 


3 


PC 2 


(0.4 ±2.4) x 10" 


-27 


(0.7 ±3.2) x 10" 27 


(0.1 ± 1.1) x 10" 


-27 


(0.5 ±1.5) x 10" 27 


3 


PC 3 


(1.7 ±4.1) x 10" 


-27 


(7.1 ±5.8) x 10" 27 


(-0.3 ± 1.7) x 10" 


-27 


(5.5 ±2.5) x 10" 27 


5 


PC 1 








(0.9 ±5.0) x 10" 


-28 


(8.87 ±0.83) x 10" 27 


5 


PC 2 








(0.1 ± 1.1) x 10" 


-27 


(0.4 ± 1.5) x 10" 27 


5 


PC 3 








(-0.1 ± 1.8) x 10" 


-27 


(5.4 ±2.7) x 10~ 27 


5 


PC 4 








(0.3 ±2.5) x 10" 


-27 


(1.9 ±2.8) x 10"" 27 


5 


PC 5 








(0.2 ±3.4) x 10" 


-27 


(3.1 ±4.1) x 10" 27 


1 


ewiMp(z) 


(0.3 ±2.2) x 10" 


-28 


(1.84 ±0.31) x 10~ 27 


(-0.1 ± 1.0) x 10" 


-28 


(1.81 ±0.16) x 10"" 27 



TABLE III: Constraints on the principal component amplitudes using simulated data for Planck and for a Cosmic Variance 
Limited experiment, assuming a fiducial energy deposition history of p ann = 1-78 x 10 _27 cm 3 /s/GeV or no energy deposition 
(pann = 0.). In the last line we also show the constraints on the ewiMp(-z) principal component presented in ijV] The 
uncertainties reported are at the 68% c.l. We show the constraints obtained when different numbers of principal components 
are varied at the same time, together with the 6 ACDM cosmological parameters. The results are reported in cm 3 /s/GeV. 



A. Constraints and forecasts 

In TablelUwe give the WMAP 7 upper limit at 95% c.l. 
obtained on the amplitude of the first principal compo- 
nent (for generic energy deposition histories) and on the 
amplitude of the principal component ewiMp(^) (the uni- 
versal WIMP energy deposition history described in jJVJ). 
For these constraints we imposed a positive fiat prior on 
the amplitude of the principal component. This physical 
assumption is convenient in order to avoid a region of pa- 
rameter space where the likelihood function is abruptly 
cut at some negative value of the principal component 
amplitude, where the recombination history calculation 
breaks down. 

Figure Q~J] shows the 1-D and 2-D contour plots from 
WMAP 7 for the first three PCs, and the forecast im- 
provements for Planck and a CVL experiment. The re- 
sults shown are obtained varying the standard cosmolog- 
ical parameters together with the amplitudes of the first 
3 principal components for each experiment. 

From these plots, it is evident that the likelihoods for 



the PC amplitudes obtained from WMAP 7 are highly 
non-Gaussian. One reason is that the WMAP 7 data 
allow sets of coefficients for the principal components 
such that the corresponding energy deposition is nega- 
tive and unphysical for some redshift. In particular, if 
the energy deposition is negative and sufficiently large in 
magnitude, this can cause the ionization history calcu- 
lation to break down. If the coefficients of the PCs are 
all large, the condition that some linear combination of 
the PCs be non-negative (or where negative, sufficiently 
small) imposes quite a non-trivial constraint on the coef- 
ficients, which is reflected in peculiar-looking boundaries 
for the favored regions. This effect is clearly visible in 
the WMAP results in Figure [TJ] for Planck the energy 
deposition history is much more constrained, and so this 
problem does not arise to nearly the same degree. 

Another issue is that adding more principal compo- 
nents to the fit does not always lead to a better recon- 
struction of the cosmological parameters and energy de- 
position history. By construction, higher PCs are less 
constrained by the data, and so the favored regions for 
their coefficients extend to much higher values, corre- 
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sponding to very large energy deposition at particular 
redshifts. This means both that the previous problem 
of unphysical energy deposition histories reasserts itself, 
and that the effects of the energy deposition become non- 
linear. 

As a consequence, the property of orthogonality be- 
tween principal components does not hold anymore, and 
unexpected degeneracies between the PC parameters can 
arise: Figure Q~J] includes a 2-D contour plot showing the 
degeneracy between PC 2 and PC 3 for WMAP 7. This 
problem affects the Planck and CVL cases only when too 
many PCs that are completely unconstrained by the data 
are included in the MCMC runs. Thus there is an opti- 
mal number of PCs to include in the reconstruction for 
each experiment: beyond this point, adding further PCs 
to the reconstruction only obfuscates the results. 

In Figure QjH we show how the bias to the cosmological 
parameters is reduced with the inclusion of more PCs, for 
both Planck and a CVL experiment. We find that the 
best results are obtained when the number of included 
PCs equals the maximum number of "measurable" PCs 
(in the sense of ^IVBp . i.e. three PCs for Planck and 
five for a CVL experiment 11 . Later PCs cannot be re- 
constructed from the data and their inclusion does not 
improve the residual bias. 

Table IIIII shows the forecast constraints (with uncer- 
tainties at 68% c.l.) on the principal components ampli- 
tudes using simulated data for Planck and for the CVL 
case, assuming a fiducial model with no energy depo- 
sition (p ann = 0) or with energy deposition described 
by constant p ann = 1.78 x 10~ 27 cm 3 /s/GcV. In these 
cases we let the amplitudes of the principal components 
assume both positive and negative values. In the ta- 
ble we report the results obtained when including dif- 
ferent numbers of principal components. The uncertain- 
ties on the amplitudes of the PCs (e.g. on PC 1) do 
not substantially change with the number of principal 
components included in the run. Table IIIII also presents 
the constraints obtained using the principal component 
ewiMp(z) described in SjV] 

Figures [TB]fT71 show the constraints on the various 
parameters from simulated Planck data with a toy- 
model energy deposition history corresponding to con- 
stant Pann(-z) = 1-78 x 10~ 27 cm 3 /s/GeV. As mentioned 
previously, while the coefficients of the early PCs are re- 
constructed well, we see that including additional PCs 
beyond the three expected to be measurable degrades 
the measurement of n s and to a lesser extend A s . 



11 Note that the number of "measurable" PCs may not reflect the 
actual number of PCs with S/N > 1 for any specific energy depo- 
sition history, which depends on the values of the a coefficients 
for that history; it is simply an estimate of the number of PCs 
that could feasibly be reconstructed from the data, from scan- 
ning over a range of models, and seems to also well describe the 
number of PCs with small enough error bars that nonlinearities 
do not cause problems. 



Figures [1811191 show the constraints on the various pa- 
rameters from a simulated CVL experiment with the 
same energy deposition history as in the Planck case. 
We now see that including at least five PCs is necessary 
to remove the bias to the cosmological parameters, espe- 
cially n s and A s but going from five to seven PCs neither 
greatly improves nor degrades the reconstruction, and is 
only significant at all in the case of n s . 



B. Application of constraints 

Once the constraints on the individual principal com- 
ponents are obtained, they can be used to set general 
bounds. Given an arbitrary energy deposition history, it 
can be decomposed into the principal component basis 
we have supplied online (see Appendix [Cj , and the co- 
efficients in that basis compared to the limits presented 
here for WMAP 7. 

As a simple example, suppose we wish to set lim- 
its on an energy deposition history with positive con- 
stant Pann, by projecting onto the constraints obtained 
for the first principal component using WMAP 7 data. 
Recall that we impose a positive prior on the ampli- 
tude of the principal component. The upper limit on 
the first principal component using WMAP 7 data is 
then 12 e x < 1.2 x l(T 26 cm 3 /s/GeV at 95% confidence. 
From Equation [9] and the discussion in £|TT1 we see that 
for a constant-pann energy deposition history we have 
£» = PannX^jLi e i( z j) : the resulting derived upper limit 
on an energy deposition history with constant positive 
Pann is Pann < 2.8 x 10~ 27 cm 3 /s/GeV at 95% confidence. 
This constraint can be compared with the one obtained 
directly by sampling the cosmological parameters with 
a positive constant p ann , which gives an upper limit of 
Pann < 2.7 x 10~ 27 cm 3 /s/GeV at 95% confidence. In a 
more general case when more than one principal compo- 
nent is considered in the analysis, the limit on p ann can 
be determined with a x 2 analysis as in Equation [TU writ- 
ing the theoretically predicted amplitude &j as a function 

Of Pann- 

For WIMP models, the WMAP 7 limit for e W iMp(z) 
in Table |TT] can be applied directly given the effective /- 
value; these were described in fjv]and are available online 
for a range of models (see Appendix [Cj) . As another con- 
sistency check, we can apply the WMAP 7 limit on the 
first principal component to the WIMP case, p an n(z) = 
eewiMP (z) ■ For WMAP 7, the coefficient ex is given by 
4.41e, giving a bound of e < 2.7 x 10~ 27 cm 3 /s/GeV. The 
true bound on e is roughly 10% stronger. 

Once Planck data are available, the same analysis can 



This value depends on the normalization of the principal com- 
ponents used. Here we employ the PCs supplied online (see Ap- 
pendix [C} , which are orthonormal when expressed as a vector 
sampled at 50 redshifts. 
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FIG. 15: Bias on the cosmological parameters using models with ACDM+ n principal components for Planck {left panel) and 
CVL {right panel) simulated data. Here, the bias is defined as (#; — 8{' td )/a0i, where 6i is the value of the mean value of the 
parameter from CosmoMC, 9{ ld is the WMAP 7 marginalized value used as fiducial for the mock data, and a{8i) is 1 — a error 
bound from the runs. Both sets of mock data assumed a constant-p an n(z) energy deposition history with p a nn = 1 X 10 _6 m 3 /s/kg 
= 1.8 x 10 -27 cm 3 /s/GeV. Note that this plot should not be directly compared to Figure [TT1 fwhich shows the bias per unit- 
normalized PC, rather than the total remaining bias after a certain number of PCs have been reconstructed from the fit). 



be redone with real data, and a very broad range of mod- 
els can then be confronted with the resulting PC-based 
limits (the only exceptions being models which have neg- 
ligible overlaps with the first few PCs and thus evade the 
constraints in this form; studying limits on such models 
will still require a separate analysis). 

Using mock Planck data, we have confirmed that the 
constraints on the first three principal components can be 
used to recover the correct limit on a particular energy 
deposition history - that is, the limit that we would ob- 
tain by directly varying the energy deposition amplitude 
for that specific model. Again taking the constant-p ann 
case as a simple example, but now using Planck mock 
data with an energy deposition history with constant 
Pann = 1-78 x 10 _27 cm 3 /s/GeV, the amplitudes for the 
first 3 principal components are: 



ei = (8.8 ± 1.5) x l(T 27 cm 3 /s/GeV (12) 
e 2 = (6.8 ± 31.5) x 10" 28 cm 3 /s/GeV, (13) 
e 3 = (7.1 ±5.8) x HT 27 cm 3 /s/GeV, (14) 



where the errors are at 68% confidence. From a y 2 anal- 
ysis of these results, the recovered energy deposition is 
given by p ann = (1.90 ± 0.50) x 10- 27 cm 3 /s/GeV at 
68% c.L. This constraint can be compared with the one 
obtained directly by sampling the cosmological param- 
eters with a constant p an n, which gives an upper limit 
of p ann = (1.90 ± 0.32) x 10- 27 cm 3 /s/GeV at 68% c.L. 
These checks confirm the validity and usefulness of the 
principal component decomposition. 



C. Bias on cosmological parameters and agreement 
with the Fisher matrix 

In Figures |2"01I2"TI we compare the error bars on the 
cosmological parameters and principal components, and 
the bias to the cosmological parameters due to assuming 
no energy deposition, for the Fisher matrix method and 
CosmoMC. We find in general that the results of the Fisher 
matrix method are in good agreement with the full like- 
lihood analysis, accurately predicting the error bars on 
the reconstructed values of the PCs and the cosmological 
parameters. For example, for the case of the Planck fit 
with three PCs, the Fisher matrix predictions for the er- 
rors on the reconstructed PCs are within 5% of the true 
values for the first two PCs, and ~ 15% different for the 
third PC, if the true energy deposition is small enough 
to lie in the linear regime. The biases to the cosmologi- 
cal parameters, when no PCs are included in the fit but 
the "true" energy deposition history is one of constant 
positive p a nn j are not quite as well matched; the Fisher 
matrix method adequately captures the directions and 
approximate sizes of the various biases, but significantly 
overpredicts the bias to A s in particular, in the example 
presented here. To obtain precise measurements of the 
biases, a CosmoMC analysis like the one we have performed 
is essential. 



VII. CONCLUSION 

Principal component analysis provides a simple and 
effective parameterization for the effect of arbitrary en- 
ergy deposition histories on anisotropies in the cosmic mi- 
crowave background. We find that for DM annihilation- 
like energy deposition histories the first principal compo- 
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FIG. 16: Constraints from simulated data for Planck on AC DAI parameters + principal components (green), ACDM+ 1 
PCs (magenta), ACDM+ 3 PCs (blue), ACDM+ 5 PCs (red). The plot shows marginalized one-dimensional distributions 
and two-dimensional 68% and 95% limits. The mock data for Planck assumed for the solid lines includes energy deposition 
with constant p ann = 1 X 10 -6 m 3 /s/kg = 1.8 x 1CP 27 cm 3 /s/GeV. The grey area shows the case of a mock data with no energy 
injection and a model AGDM+ principal components. Only the cosmological parameters are shown. 



nent, describing the bulk of the effect, is peaked around 
z ~ 500 — 600, at somewhat lower redshift than previ- 
ously expected; the later principal components provide 
corrections to this basic weighting function. 

The principal components, derived from a Fisher ma- 
trix approach, are stable against a wide variety of per- 
turbations to the analysis, including choice of code cal- 
culating the ionization history, additions to the usual set 
of cosmological parameters, the inclusion or exclusion of 
ionization on helium, the range of included multipoles, 
and the choice of binning. (For further discussion, see 
Appendix [5]) The one significant potential change to 
the PCs arises from how deposited energy is attributed 
to additional Lyman-a photons: we have showed the ef- 
fect of on one hand neglecting this channel, and on the 
other of assuming that all the energy attributed to "exci- 
tations" is converted into Lyman-a, which should bracket 
the true result. We eagerly await a more careful analysis 



of this problem. 

Within the Fisher matrix formalism, it is straightfor- 
ward to take into account degeneracies with the standard 
cosmological parameters. We have presented predictions 
for the (significant!) biases that would arise in Planck as 
a result of falsely assuming energy deposition to be zero, 
for each of the principal components. We have confirmed 
the previously noted degeneracy between energy deposi- 
tion and n s , and to a lesser degree with A s , uj^ and u c , 
in WMAP data; since our analysis decomposes the bi- 
ases according to the principal components that generate 
them, it is now trivial to compute the biases to the cos- 
mological parameters for any arbitrary energy deposition 
history, in WMAP or in mock Planck data. 

For a wide range of energy deposition histories, span- 
ning models of dark matter annihilation and decaying 
species where annihilation or decay can begin or end 
abruptly on characteristic timescales shorter than the age 
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FIG. 17: Constraints from simulated data for Planck on ACDM+ 1 PCs (magenta), ACDM+ 3 PCs (blue), ACDM+ 5 PCs 
(red). The plot shows marginalized one-dimensional distributions and two-dimensional 68% and 95% limits. The mock data for 



Planck assumed for the solid lines includes energy deposition with constant p a 



with the PC coefficients shown in units of m /s/kg. Only the principal components and n a are shown. 



1 x 10- b m7s/kg = 1.8 x 1(T 27 cm7s/GeV, 



of the universe, the coefficients of up to three principal 
components are potentially measurable by Planck, for en- 
ergy deposition histories satisfying 95% confidence limits 
from WMAP, opening up the exciting possibility of dis- 
tinguishing different models of energy deposition. For a 
CVL experiment, up to five coefficients are measurable. 

For the "standard" WIMP annihilation case, principal 
component analysis on a large set of WIMP models yields 
a single principal component ewiMp(z) that describes the 
effect on the GYs of all the models very well; any model is 
then parameterized simply by the coefficient of ewiMp(z) 
(or equivalcntly, effective /). Our analysis confirms pre- 
vious statements in the literature, and we have provided 
this "universal f(z) n curve for future WIMP annihilation 
studies. 

We performed an accurate MCMC analysis of current 
WMAP 7 data to impose constraints on the measurable 
principal component amplitudes, and to forecast con- 
straints for future experiments such as Planck or a CVL 



experiment. We find good agreement with the Fisher ma- 
trix analysis, although the MCMC analysis is required to 
accurately predict the biases on the cosmological param- 
eters. We have illustrated how it is possible to recover 
the constraints on an arbitrary energy deposition history 
from the constraints on the amplitudes of the principal 
components. The reconstructed constraints are in very 
good agreement with the constraints obtained by directly 
sampling a specific energy deposition history, confirming 
the validity and usefulness of the principal component 
decomposition. 
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FIG. 18: Constraints from simulated data for a cosmic variance limited experiment on ACDM+ principal components 
(green), AC DM + 1 PCs (magenta), ACDM+ 3 PCs (blue), ACDM+ 5 PCs (red), ACDM+ 7 PCs (black). The plot shows 
marginalized one-dimensional distributions and two-dimensional 68% and 95% limits. The mock data for the CVL experiment 
assumed for the solid lines includes an energy deposition history with constant p ann = 1 x 10 _6 m 3 /s/kg = 1.8xl0 -27 cm 3 /s/GeV. 



The grey area shows the case of a mock data with no energy injection and a model ACDM- 
the cosmological parameters are shown. 



principal components. Only 



DE-FG02-90ER40542. This research made use of the Astronomy User's Library at Goddard 13 . 
NASA Astrophysics Data System (ADS) and the IDL 

13 Available at http://idlastro.gsfc.nasa.gov 
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FIG. 19: Constraints from simulated data for a cosmic variance limited experiment on ACDM+ 1 PCs (magenta), ACDM+ 
3 PCs (blue), ACDM+ 5 PCs (red), KCDM+ 7 PCs (black). The plot shows marginalized one-dimensional distributions and 
two-dimensional 68% and 95% limits. The mock data for the CVL experiment assumed for the solid lines includes an energy 
deposition history with constant p an n = 1 X 10 _6 m 3 /s/kg = 1.8 x 10~ 27 cm 3 /s/GeV, with the PC coefficients shown in units 
of m 3 /s/kg. Only the principal components and n s are shown. 
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FIG. 20: Error bars for the cosmological parameters and coefficients of the principal components, in mock Planck data, 
simulated assuming a constant-p ann (z) energy deposition history. The points marked "large p ann " have pann — f x 10 _6 m 3 /s/kg 
= f.8 x 10 -27 cm 3 /s/GeV, while for the points marked "small p ann " the value is a factor of 10 lower. In the left panel the 
fit to the data is performed using only the standard six cosmological parameters, whereas in the center panel and right panel 
three PCs are also included in the fit. For ease of comparison, the green triangles in the center panel reproduce the black 
asterisks in the left panel. For the PCs, the dimensionless uncertainties plotted here should be multiplied by p arm to obtain the 
uncertainties on the e« coefficients. The uncertainties on the cosmological parameters have been divided by the central values 
of the respective parameters. 
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FIG. 21: Biases to the cosmological parameters in mock 
Planck data, simulated assuming a constant-p ann (,z) energy 
deposition history with p ann = 1 x 10~ 6 m 3 /s/kg = 1.8 x 10~ 27 
cm 3 /s/GeV. The fit to the data is performed using only the 
standard six cosmological parameters. 
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FIG. 22: An example of the dependence of the principal components on the number of bins and the width of the Gaussians, 
for the first, third and sixth principal components in Planck, after marginalization over the cosmological parameters (the later 
principal components are shown as they are more sensitive to small changes in the binning). Black lines indicate the case with 
100 linearly spaced bins and o — bin width/2; in the left (right) panel, red and blue dots indicate respectively the case with 15 
(50) linearly spaced bins and a = bin width/2, and the case with 15 (50) linearly spaced bins and a = bin width/4. 

Appendix A: Validation of the PCA method 

In this appendix we discuss a number of specific issues that might affect the results of the principal component 
analysis, and motivate our default choices. 



1. Choice of binning 



As discussed previously, to approximate r5-functions in energy deposition, as a default we employ (for the 
annihilation-like case) 50 linearly-spaced redshift bins covering the redshift range from z — 80 to z = 1300. Fig- 
ure [HI shows the effect of changing the Gaussian width and the number of bins for some sample PCs; doubling the 
number of bins or changing the Gaussian width by a factor of 2 has no effect on the PCs, although reducing the 
number of bins to 15 does affect the PCs slightly, especially at the lowest and highest redshifts. 

In f|TT]we briefly discussed the choice of log vs linear binning, preferring to use the former for decay-like scenarios 
and the latter for annihilation-like scenarios. In the annihilation-like case, the choice of linear-spaced bins in redshift 
was not inevitable; log-spaced bins seem equally natural. The choice of log or linear binning is somewhat subtle, as 
it affects whether or not two energy deposition histories are considered orthogonal, and how much different redshifts 
contribute to the norm of a particular energy deposition profile. Consequently, the two choices give rise to different 
sets of principal components and eigenvalues, and a "generic" energy deposition history is somewhat different between 
the two cases. However, the eigenvalues of the first several PCs are quite similar, and using log binning instead of 
linear binning does not significantly affect the results of EIIVBI in both cases a similar (small) number of PCs are 
sufficient to describe a very broad class of energy deposition histories, for the purpose of future experiments. For 
individual models, of course, the number of measurable PCs depends on the choice of basis: a p(z) curve which 
happens to be very well described as a linear combination of the first two PCs in the log-binned case may not be 
nearly so well described by the first two PCs in the linear-binned case, or vice versa. 

In Figure [3] we showed the Planck PCs for both linear and log binning in the annihilation case. The first PC is 
peaked at higher redshift in the log-binned case: this is to be expected, since in the log-binned case the amount of 
energy deposited per bin for a constant p(z) has an extra (1 + z) scaling relative to the linear case (simply due to 
the wider bins at high redshift). It can be shown that this simple rescaling largely (but not completely) describes the 
difference between log and linear binning for the first PC, but the later PCs shift by a larger amount and in more 
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FIG. 23: The dependence of a selection of PCs on the range of i included in the analysis, for WMAP 7 (left) and Planck (right). 
complex ways, as each must be orthogonal to all previous PCs. 



2. The effect of noise and maximum I 

As a default, we include I — 2. .2500 in our analysis, for both temperature and polarization. Ps above 1000 do not 
noticeably affect the principal components for WMAP, but have significant effects for Planck (and the CVL case), 
especially in the higher PCs. Figure |2"51 shows the effect on a selection of PCs on changing ^ max . The results also 
depend on the estimated sensitivity of the experiment. 



3. Choice of ionization history calculator and recombination corrections 



Recent improvements in the detailed treatment of recombination have motivated new codes to compute the ioniza- 
tion history, in particular CosmoRec [l5| and HyRec [lfjl ]. The inclusion of these additional recombination corrections 
can modify the effect of energy deposition on the ionization history and the CMB, thus shifting the derivatives and 
modifying the principal components. However, comparing the principal components obtained using RECFAST 1.5 to 
those obtained using CosmoRec, we hnd that there is essentially no difference. The overall normalization of the effect, 
for a given energy deposition history, is slightly smaller (by a few percent) in CosmoRec as compared to RECFAST 1.5, 
but as shown in Figure [24] the shapes of the PCs are unaffected. This shift in normalization will not change the 
detectability of the maximum WMAP -allowed signal in Planck, and will only very slightly modify the amount of 
energy deposition required to produce such a signal. 



4. Helium and Lyman-a 



Ionization of helium also contributes to the population of free electrons, and is included automatically in CosmoRec: 
however, some previous studies have neglected it, focusing only on the H contribution. We include helium by default 
throughout, but show in Figure [24] the (negligible) effect on the PCs of leaving it out. Similarly to the choice of 
ionization history calculator, including helium increases the total effect of a given energy deposition history, but does 
not significantly modify the redshift dependence and hence the PCs. 

A more difficult question is the effect of additional Lyman-a photons. The energy from dark matter annihilation that 
is deposited to the gas is partitioned between excitation, ionization and heating. In the analysis of Q , which we have 
followed, only the latter two processes are included; excitations are assumed to have no significant effect. However, 
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FIG. 24: The first three principal components for Planck, after marginalization, computed using RECFAST 1.5 and CosmoRec. 
In the baseline case (as in CosmoRec), ionization of helium is included but injection of Lyman-a photons is not. We also show 
the effects of including a contribution to Lyman-a photons, and neglecting helium. The effect of helium ionization on the PCs 
is negligible because it is approximately a redshift-independent effect. 



it was pointed out in [36jj that the increased population of Lyman-a photons from excitations - or equivalently, the 
higher fraction of H in an excited state - facilitates ionization and thus indirectly increases the ionization fraction. 

We studied the effect of including Lyman-a photons, following |3g, HH and assuming all the energy partitioned 
into excitation produces Lyman-a photons. The true picture is probably more complicated, but its study lies beyond 
the scope of this work, and this choice and our baseline case with no additional Lyman-a photons should bracket 
the true solution. Figure [24] shows the effect on the PCs of including the Lyman-a contribution. The corresponding 
effect on the S/N for the various PCs is indicated by the hatched region in Figure [TU1 the bound on constant p ann is 
strengthened by ~ 7 — 10% (depending on the experiment). 



5. Additional cosmological parameters 

Since the primary effect of dark matter annihilation is an ^-dependent damping of the temperature anisotropics, 
one might ask if it is degenerate with other cosmological parameters that have a similar effect: for example, including 
running of the scalar spectral index could better mimic the profile of the damping with respect to increasing 
the number of massless neutrino species N e g also suppresses small-scale anisotropy, and the primordial helium mass 
fraction Yp has degeneracies with N e g (56[. Accordingly, we have redone the principal component analysis including 
each of these parameters separately, as well as all of them simultaneously. 

In all cases, we find only negligible shifts to the principal components. The forecast constraint on constant p ann 
from Planck is weakened by ~ 7% if running of the spectral index is included, and by ~ 12% if all three parameters 
are used. For WMAP 7, the Fisher-matrix-estimated limit weakens by ~ 11% and ~ 13% in these two cases. 



Appendix B: Review of marginalization and biases 



1. Marginalization as projection 

Marginalization over some parameters in the Fisher matrix is equivalent to a projection of the vector space spanned 
by all the parameter perturbations. Another way of saying this is that given the transfer matrix (Equation [2]) for 
the parameters of interest, then the marginalized Fisher matrix for those parameters is formed by contracting the 
projected transfer matrix with the covariance matrix £ . 

We illustrate these statements for the situation considered in this paper. The total vector space consists of the 
SCis spanned by energy deposition and cosmological parameter perturbations, and we wish to marginalize over the 
cosmological parameters. Recall that the full Fisher matrix is 



( F e F v \_f TjE- 1 ^ TjE _1 T c \ 
{F? F c )-\ Tf E-iT e 7?£-iT c ) 
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where T e and T c are transfer matrices mapping energy deposition histories and cosmological parameter perturbations, 
respectively, to the space of SC/s. For ease of notation we suppress indices. Using Equation [BJ the marginalized 
Fisher matrix can be written as 

F = F e — F V F- X F% = (PT e f ST 1 (PT e ) (B2) 

where P = (l — T c (F c ) 1 T^S -1 ^ satisfies P 2 = P and so is a projection operator. P projects out any component 

of SCe which can be effectively absorbed by a change in the cosmological parameters. This can also be seen if we act 
with P on a generic perturbation in the cosmological parameters: PT c S9 a = 0. Accordingly, any perturbation to the 
fiducial CMB model can be written as 

Sd = (1 - P)Sd + PSd = SC l e l + 8Cj. (B3) 

where _L and 1 1 mean perpendicular or parallel to the cosmological parameter perturbations. 

Thus PT e is a projected transfer matrix, taking energy deposition histories to the subspace of SC^s which are 
orthogonal to cosmological parameter perturbations. This projection depends on the noise parameters in S -1 , since 
the notion of orthogonality depends on our definition of norm. Intuitively, if an energy deposition and cosmological 
parameter perturbation have very similar effects at low I but are different at high £, then the projection operator may 
give a very small 8Cj~ in the case of WMAP and a larger 8Cj- in the case of Planck. 

The eigenvectors ej of F (or principal components) with the largest eigenvalues are correspondingly those with the 
largest measurable 8C^. The ei also map to an orthogonal vector space of SCjh's. The unprojected hi — T e ti are in 
general not orthogonal because hi ■ hj — (T e ei) T E _1 (T e ej) = ej F e tj. However, defining hj~ = PT e ei, we have 

hi ■ hj = (PT e e 4 ) T I]- 1 (Fr e e J ) = ejF ej = ^A 4 . (B4) 

Because the hi~ are orthogonal, we can extract the coefficients of the principal components from a generic Ri G {SCe}: 

£l = UpT^f^'Ri = (B5) 



2. Biases to cosmological parameters 

The parallel components from energy deposition, £C] = (1 — P)T e eiEi, correspond to the biases to the cosmological 
parameters. Suppose there is some energy deposition (si ^ 0), but it is not included in the fit to the data. Then the 

best fit cosmological parameters will absorb any SC e from the energy deposition, and the measurements of the true 

cosmological parameters would be biased by the amount it takes to produce 8C\ : 

SC l e l = T c S6 a - (1 - P)T e e l£l . (B6) 
Multiplying both sides by (Fj-^jE" 1 , we have 

S6 a = F^F^em- (B7) 



Appendix C: Web files description 

To facilitate the analysis described above, we provide the PC vectors defined above on our web page 14 . The files 
are provided in two formats: as a the Flexible Image Transport System 15 (FITS) binary table, and as ASCII plain 
text files. The file formats are described briefly in this appendix, and in more detail on the web page. 



http: / /nebel. rc.fas.harvard.edu/epsilon 
http:/ /fits. gsfc.nasa.gov 
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FITS files 



Each FITS file has the following format: 



Array [2500] 
Array [50] 
Array [50, 5] 
2.0000000e-27 
Array [5] 

Array [2500, 3, 5] 
PC_POWSPEC_PERP Array [2500, 3, 5] 



L 

REDSHIFT 

PC_EIGENVECTORS 

EPSILON 

EIGENVALUES 

PC_POWSPEC 



- multipole index (2.. 2501) 

- redshift (z) 

- PC as a function of redshift 

- energy deposition, see Section IIIB of paper 

- PC eigenvalues 

- PCs projected into CI space, TT, EE, TE 

- same, but projected onto the space perpendicular to the 
cosmological parameters 



We also provide the universal ewiMP curve derived in ifVl following the same notation, and a list of efficiency 
coefficients for the WIMP models studied in 11011 . 



2. ASCII files 

We also provide the mappings of the PCs into <5C/s in ASCII files, one per principal component. File names are 
e.g. 

epsilon_ann_Planck_PC01 . dat 
epsilon_ann_Planck_PC01_perp . dat 

etc., with one file for each PC, choice of binning, and experiment (CVL, Planck and WMAP 7). Each of these 
ASCII files contains 4 columns: £,PC T t,PCee,PC T e- 

For each experiment and choice of binning, one additional ASCII file contains the ordered eigenvalues for the 
principal components, and another holds the ordered eigenvectors / PCs as functions of redshift. These files are 
named e.g. 

epsilon_ann_Planck_PC_eigenvalues . dat 
epsilon_ann_Planck_PC_eigenvectors . dat 



3. Units 



The PC vectors PC_P0WSPEC and PC_P0WSPEC_PERP are the changes in £(£ + l)Ci/2i: (before and after projecting 
out the cosmological parameters, respectively) corresponding to an energy deposition history given by £i = e for all 
i (where the £i are dimensionful coefficients defined in Equation [9]). Our convention, described in £|III B[ is that the 
Ci have units of /zK 2 , the principal components are dimensionless, energy deposition from annihilations has units 
cm 3 /s/GeV, and energy deposited by decays has units of s _1 (see Equation [1} . The energy deposition parameter e 
(labeled EPSILON in the files) is fixed at 2 x 10~ 27 cm 3 /s/GeV for the annihilation case, and 1 x 1CT 26 s _1 for the 
decay case. 
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